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Fluctuations of the electromagnetic field produced by quantized matter in external electric field 
are investigated. A general expression for the power spectrum of fluctuations is derived within the 
long-range expansion. It is found that in the whole measured frequency band, the power spectrum 
of fluctuations exhibits an inverse frequency dependence. A general argument is given showing that 
for all practically relevant values of the electric field, the power spectrum of induced fluctuations 
is proportional to the field strength squared. As an illustration, the power spectrum is calculated 
' explicitly using the kinetic model with the relaxation-type collision term. Finally, it is shown that 

the magnitude of fluctuations produced by a sample generally has a Gaussian distribution around 
(***«) , its mean value, and its dependence on the sample geometry is determined. In particular, it is 

■ demonstrated that for geometrically similar samples, the power spectrum is inversely proportional 

^ ' to the sample volume. Application of the obtained results to the problem of flicker noise is discussed. 
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, I. INTRODUCTION 



As is well-known, power spectra of voltage fluctuations in all conducting materials exhibit a universal profile in 
the low-frequency limit, which is close to inverse frequency dependence. Fluctuations characterized by the power 
spectrum of this type are called usually 1//, or flicker, noise. Although this noise dominates only at low frequencies, 
experiments show the presence of the l//-component in the whole measured frequency band up to 10 6 Hz. Despite 
""^ \ numerous attempts, no lower frequency bound for the 1// law has been found. In addition to that, it is generally 
accepted that l//-noise produced by a sample is universally characterized by the following properties 1) it is (roughly) 
inversely proportional to the sample volume, 2) it is Gaussian, and 3) its part induced by external electric field is 
proportional to the field strength squared. 

A number of mechanisms has been suggested to explain the origin of 1/ /-noise lj . There is a widespread opinion 
that this noise arises from resistance fluctuations, which is quite natural taking into account the property 3) mentioned 
CJ^ 1 above. Indeed, for a given current through the sample, the mean squares of voltage and resistance fluctuations are 
\ proportional to each other, the current squared being the proportionality coefficient. It has been proposed that 
the resistance fluctuations possessing the other properties of flicker noise might result from temperature fluctuations 
0, fluctuations in the carrier mobility 0, 3 or in the number of carriers caused by surface traps @. All these 
models, however, have restricted validity, because they involve one or another assumption specific to the problem 
under consideration. For instance, assuming that the resistance fluctuations arise from the temperature fluctuations, 
one has to choose an appropriate spatial correlation of these fluctuations in order to obtain the desired profile of the 
power spectrum. Similarly, the model of Ref. requires specific distribution of trapping times. In addition to that, 
the models proposed so far reproduce the 1 / /-profile only in a restricted range of frequencies, require an appropriate 
normalization of the power spectrum, etc. At the same time, ubiquity of flicker noise and universality of its properties 
suggest existence of a simple and universal, and therefore, fundamental origin. It is natural to look for this reason 
in the quantum properties of charge carriers. In this direction, the problem has been extensively investigated by 
Handel and co-workers Handel's approach is based on the theory of infrared radiative corrections in quantum 
electrodynamics. Handel showed that the 1// power spectrum of photons emitted in any scattering process can be 
derived from the well-know property of bremsstrahlung, namely, from the infrared divergence of the cross-section 
considered as a function of the energy loss. Thus, this theory treats the l//-noise as a relativistic effect (in fact, the 
noise level in this theory ~ a(Av) 2 / c 2 , where a is the fine structure constant, Av velocity change of the particle being 
scattered, and c the speed of light). It should be mentioned, however, that the Handel's theory has been severely 
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criticized in many respects 7, 8J. 

In Refs. P. Il0| . the role of quantum effects is considered from a purely nonrelativistic point of view. In Ref. [jj, 
quantum fluctuations of the electromagnetic field produced by elementary particles are investigated, and it is shown, 
in particular, that the correlation function of the fluctuations exhibits an inverse frequency dependence in the low- 
frequency limit. This result was applied in Ref. [Tpj to the calculation of the power spectrum of electromagnetic 
fluctuations produced by a sample. It was proved, in particular, that the power spectrum possesses the properties 1), 
2) of flicker noise, mentioned above. As to the property 3), it was argued in Ref. |10| that this requirement is also 
met. The argument was based on the assumption of analyticity of the electron density matrix with respect to the 
external electric field. However, this assumption is not valid in general. Thus, the issue concerning the influence of 
external field is left open. 

The purpose of the present paper is to investigate the influence of the external electric field on quantum electro- 
magnetic fluctuations in detail. We will show that for all practically relevant values of the field strength, the power 
spectrum of induced fluctuations is proportional to the field strength squared indeed. 

Inclusion of an external field lowers the system symmetry, therefore, our first problem below will be to generalize 
the results obtained in |lf)j | for spherically-symmetric systems to systems with axial symmetry. This is done in Sec. Mil 
First of all, we prove in Sec. MI XI that the low- frequency asymptotic of the connected part of correlation function is 
logarithmic. This result was proved already in [T(il |. Although the proof does not rely on the system symmetries, we 
give an independent and more simple and accurate proof of this important fact. Contribution of the disconnected 
part is calculated in Sec. MI Bl and is found to exhibits an inverse frequency dependence, thus dominating in the 
low- frequency limit. Because of the lower system symmetry, this calculation is much more complicated than in the 
case considered in |Io|. The obtained expression for the power spectrum is analyzed in Sec. II VI where a general 
argument is given showing that the field-induced noise is quadratic in the field strength, which is then illustrated 
using the simplest kinetic model with the relaxation-type collision term. SectionlVlsummarizes the results of the work 
and states the conclusion. 



II. PRELIMINARIES 



Let us consider electromagnetic field produced by a classical resting particle with mass m and electric charge e. It 
is described by the Coulomb potential 

Ao = ^, A = 0. (1) 

In quantum theory, this form of the electromagnetic potential is reproduced by the mean fields (in|Ao|in), (m|A|m) 
calculated far away from the region of particle localization. If the 3-vector of the mean particle position is denoted 
by x , and that of the point of observation by x, then the latter condition means that \x — x \ ^> D, where D 
is a characteristic length of the particle wave packet spatial spreading (for instance, the variance of the particle 
coordinates). As a result of the quantum evolution according to the Schrodinger equation, D increases in time, 
thus leading to a dispersion of the electromagnetic field produced by the particle. Furthermore, because of the 
quantum indeterminacy in the particle position, the field fluctuates. The correlation function of the fluctuations is 
conventionally defined by 

C^{x,x') = -{in\A^{x)A v {x') + A v {x')A^{x)\ixi) - (in|i M (ar)|in) (in|i„(a;')|m) , (2) 

where x and x' are the spacetime coordinates of two observation points. Of course, this function is dispersed, 
too. Our aim below will be to investigate low-frequency properties of this dispersion. It is clear that the condition 
|cc — a? 1 3> D is irrelevant in this investigation, because the low- frequency asymptotic of the power spectrum of 
correlations is determined largely by the late-time behavior of the function C fiu (x,x'), where D is unbounded (for a 
free particle state, and for large times t, D is a linear function of t). 

As in Ref. [l(J, we will work within the long-range expansion of the correlation function, which is a convenient 
tool for extracting the leading term of the correlation function. Let us briefly recall the reasons justifying application 
of this expansion. The function C^ u (x,x') can be represented as a power series in the ratios l c /R and D/R, where 
l c = h/mc is the Compton length, and R is either r = \x — xq\ or r' = \x' — Xq\. However, as mentioned above, the 
ratio D/R cannot be considered small as long as one is concerned with the low- frequency behavior of correlations. 
We overcome this problem by going over to the momentum space, and work with an expansion in powers of l c and 
pj D, where p is the 3-momentum transfer to the particle, and D = J (q 2 ) is the variance of the particle momentum. 
Unlike D, the quantity D is time independent (for free particle states), so the expansion is valid for all times. By the 
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order of magnitude, the relevant values of the momentum transfer \p\ ~ H/r, and therefore, validity of the expansion 
in powers olp/D requires only that 

rD»/l, r'D^>h. 

All subsequent considerations are carried out under these conditions. 

We recall also that, as was shown in Ref. Q, the leading term of the correlation function is of zcroth order in 
l c , hence, in the units h = c = 1 used from now on, it can be identified as the limit of the correlation function for 
m — > oo. It should be emphasized that this identification is only formal, in particular, it does not mean that the 
results obtained below apply only to heavy particles. 

According to Eq. J2J), in order to find the correlation function of the electromagnetic fluctuations, one has to calculate 
the in-in expectation values of the field operators as well as of their products. In Refs. |9lll0|. this was done using 
the Schwinger-Keldysh formalism. In particular, it was proved in Ref. 10] that the leading low-frequency term of 
the correlation function is contained entirely in its disconnected part [the second term in Eq. (J2J]. The connected 
part of the correlation function was taken in Ref. 10] in the nonsymmetric form (ui\A ti (x)A l/ (x')\'ui). Although the 
proof given there can be carried over to the present case, we will give an independent and more accurate proof of this 
important fact, which avoids complications of the Schwinger-Keldysh method. 

Note first of all, that for all values of t, t', the connected part of the correlation function can be rewritten as 

^{m\A^x)A v {x') + A v {x')A^x)\m) = ^n\T{A^x)A v (x')} + f{A^x)A v (x')}\\n) , 

where the operation of time ordering T (T) arranges the factors so that the time arguments decrease (increase) from 
left to right. Furthermore, for a one-particle state under stationary external conditions, the state vector |in) can be 
substituted by the vector |out), up to a phase factor. In the tree approximation, this factor is equal to unity, therefore, 
one can write, taking into account that A^ is Hermitian, 

^m\A^x)A u {x') + A u {x')A^{x)\in) = i(out|T{i M (x)i,,(x')}|m) 

+ ^(m\f{A^x)A v (x')}\out) = Be{aat\T{A^(x)A v (a/)}\m) . (3) 

The latter quantity can be calculated by applying the usual Feynman rules. 

Let us assume, for simplicity, that the field-producing particle is described by a complex scalar <f>. It is not difficult 
to show actually that the results derived below are valid for particles of any spin. This is because in the long-range 
limit, the value of the electromagnetic current is fixed by the standard normalization conditions for the one-particle 
state, which are universal for all particle species. Let the gradient invariance be fixed by the Lorentz condition 1 

G = d"A^ = . (4) 

Then the action of the system takes the form 

s[$\ = s m + s gi [$>} , 

So[$] = Jd 4 x{(d^* +ieA^*){d^-ieA^)-m 2 ^<i ) }~- A Jd i xF^F^ , 

Sgf[*] - -\Jd A xG\ F llv =d ll A v -d v A ll , (5) 

where the Feynman weighting of the gauge condition is assumed. The tree diagrams generated by this action, which 
contribute to the right hand side of Eq. (j2J, are depicted in Figs. El 

To complete this section, let us define the power spectrum function of fluctuations. We are concerned with corre- 
lations in the values of the electromagnetic fields measured at two distinct time instants (spatial separation between 
the observation points, \x — x'\, is also kept arbitrary). Accordingly, fixing one of the time arguments, say, t', we 
define the power spectrum function as the Fourier transform of C^ u {x, x') with respect to (t — t') 

+oo 

C lw (x,x / ,t',uj)= f dtC^ x'y^-^ . (6) 



1 The proof of the gauge-independence of the leading contribution, given in , does not rely on the symmetry properties of the particle 
wave function, and hence carries over to the present case. 



4 



III. EVALUATION OF THE LEADING CONTRIBUTION 



Evaluation of the low- frequency asymptotic of the correlation function proceeds in two steps. First, we will prove 
in Sec. IIII Al that the low-frequency asymptotic of the connected part of C M „ [the first term in Eq. is logarithmic. 
The contribution of the disconnected part will be calculated in Sec. IIII Bl It will be shown that this contribution 
exhibits an inverse frequency dependence, and thus dominates in the low-frequency limit. 



A. Low-frequency asymptotic of the connected part of correlation function 

Before going into detailed calculations, let us first exclude the diagrams in Fig. [2J which do not contain the H° 
contribution. It is not difficult to see that these are the diagrams without internal matter lines, i.e. P^c) in the present 
case. Indeed, this diagram is proportional to the integral 

p ik(x — x') 

d A k-, 



k 2 (k — p) 2 

which does not involve the particle mass at all. Taking into account that each external matter line gives rise to the 
factor (2e q ) _1 / 2 , where e q = \J m 2 + q 2 « m, we see that the contribution of diagram[2Ic) is proportional to 1/m. 
Hence, on dimensional grounds, this diagram is proportional to h. 
The contribution of diagrams GIa),GIk) has the form 



I„ v (x,x') = ie 2 [[ d i zd i z'D a (x,z) 



^(2)^£>(2,z')^<M2') 



D u (z',x'), (7) 



where 



d 4 k e -ik(x-y) 

D(x,y) 



(2tt) 4 m 2 - k 2 - iO ' 
D° ee D| m=0) (8) 



and (f)Q is the given particle state. Introducing the Fourier transform of I tJlV (x, 

L v {x,x',t',u)= f dtl„ v {x, x'y"^ , 



and going over to the momentum space, one finds 



= [f *« »'%>°fa + ^ e-^'-^'Up, g, x - u,) , (9) 

J J [2iry (27r) J y?2e q 2e q+p 



where 



/d 3 k i 
__^k( x - x ) { 2 q ^ + k ^(2q u + k v + p v ) 

x D°(k)D(q + k)D°(k-p)\ k0 ^ . (10) 

Here is the particle 4- momentum, and a(q) its momentum wave function at some time instant to. The function 
a(q) is normalized by 

' rq ' a ( g )| 2 = l, (11) 



(2tt) 3 
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and is generally of the form 

a(q) = b(q)e- i ^ , (12) 

where Xq is the particle mean position, and b(q) describes the momentum space profile of the particle wave packet. 

Let us now show that the low-frequency asymptotic of J^ v is logarithmic. We note, first of all, that in the long- 
range limit, the 4-momenta ku,,Pu m the vertex factors can be neglected in comparison with q^, because the leading 
contribution comes from integration over small k^. By the same reason, the factor exp{ife(a; — x')} can be set equal 
to unity. Next, introducing the Schwinger parameterization of the propagators, we rewrite J M „ as 

f d 3 k 

Jp V (p, q,u) = Aq^ J ^—3 exp {i (x[uj 2 - k 2 } + y[uj 2 - k 2 + 2mto - 2(kq)\ 

+* - ^ - <* - -m - ^//jf *** (^f - {-^} 

x exp {i (xui 2 + y[uj 2 + 2mu] + z[[u - p ) 2 - p 2 ]) } . 
Changing the integration variables y — > xy, z — > xz, and integrating over x gives 



J,Mu{p,q,u) = dydz {{yq - zpf 



0\2 „2l\ 1 --V 2 



+ (1 + y + z) (lu 2 + y[uj 2 + 2muj] + z[(uj - p ) 2 - p 2 ]) } 

The singularity of the latter integral at lu — comes from integration over small y, z, therefore, it is the same as the 
singularity of the integral 

= ^JJ dydz{y 2 q 2 ~2yz(qp)+LO 2 + y[LO 2 + 2mLo} + z[(u;-p ) 2 -p 2 }y 3/2 

%<lv f°° d V r 2 2 , 2 , r 2 , o U^ 1 / 2 

2tt J (u! - p ) 2 p 2 - 2y(qp) 

= / 7 0^2 2 — o — 1 — \\y 1 +!+# + 2m ) 

2tt J {uj-p u Y p 2 2ujy{qp) 

After the transformations performed, the singularity of the last integral for lu — > reappears at y — > 00, and hence, 
it coincides with the singularity of the integral 

in 1 \ %iv f°° d v 1 



27r|q| J y (u - p ) 2 - p 2 - 2ujy{qp) 

ln y 



q ^ v In-. (13) 



2ir\q\p 2 w 



27r|q|{(tj — p ) 2 — p 2 } (lu — p ) 2 — p 2 — 2ujy(qp) 

It is not difficult to verify that the obtained expression agrees with the results of Sec. Ill A of Ref. ^j- Thus, we 
have proved that the connected part of the correlation function diverges for lu — > only logarithmically. 

B. Low- frequency asymptotic of the disconnected part of correlation function 

Let us turn to the disconnected part of the correlation function. To find its Fourier transform, we have to evaluate 
the integral 



+ 00 



^(inlA^aOHe^-*') = ^-^'/(r, w) , (14) 



where 

+00 



r = X — Xq ■ 
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To the leading order of the long-range expansion, b(q + p) can be substituted here by 6(g) . Using also 

2m 2m 

gives 

The function 6(g) will be assumed to possess the symmetry of the external field. Thus, in the presence of a homo- 
geneous external electric field, the function 6(g) is axially-symmetric; taking z axis in the direction of the field, one 
has 6(g) = b(q±,q z ), where q±^ = \J q 2 — q 2 is the transverse component of the particle momentum. In this case, I is 
a function of £ = r z and p — yV 2 — £ 2 , j(r, w) = I(p, C, w), and hence, averaging over transverse directions, one can 
write 

where is the angle between p and z axis, and Jo is the Bessel function. To evaluate the p integral it is convenient 
to introduce a spherical coordinate system, with the polar axis pointing in the direction of the vector q. Let the 
azimuthal and polar angles of p in this system be denoted by <p, 9, and those of z axis by <p' ,6' , respectively. Then 



and 



cos 4> = sin 9 sin 9' cos(</? — cp') + cos 9 cos 9' , 



%C,w = 7^ / 7^4 /3 fo]., 9, 2 / d0sin0 / rfy ° f g ==- , 15 
(2tt)^ J (27r)' 3 7 7 V? cos 2 6» + 2mw 



u = — I q\ cos + \J q 2 cos 2 + 2mw , 

where it is assumed that w > 0. Since the integrand in this formula depends on the difference <p — <p', the result of 
integration over ip is independent of ip' . Setting the latter equal to zero, and taking <f> as the integration variable yields 

2tt 6+6' 

ApUup^)***** = 2 / i<t> - ^MvptinW™* . (16) 

J ^[cos - cos(0 + 0')][cos(0 - 9') - cos 0] 

It is seen from Eq. (|15|l that the singularity at u> — comes from integration over 9 w 7r/2. Without changing the 
singular contribution, therefore, one can set = 7r/2 in the integral l|16|l . which after simple transformations takes 
the form 

r Jo (upy/l — sin 2 9'v 2 ) cos (uv(s'm9') 
3 = 4/ dv^ ? =J= . (17) 





In the long-range limit, the argument of the Bessel function is large, so one can use the asymptotic formula 

T 



■h{z) = \ — cos (z - y 

V 7TZ V 4 



Thus, the above integral involves rapidly oscillating trigonometric functions, and the leading contribution comes from 
integration around the point of stationary phase of the integrand. Assuming £ > 0, and decomposing the product of 
cosines into a sum, one sees that the range of integration contains one such point, if C < rsin#', namely 

c 

vq = — ^7 , 
r sin o 
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and hence 



[n \ exp li (up\/ 1 — sin 2 9'v 2 + uv( sin 9' )} 

:i = 2d- Re / dv L V = 

Y (1 — v 2 )up\f 1 — sin 2 9'v 2 

( . ( UT 3 . 2 , 2 7T 

I exp < i iur — —2 sin 9 [v — vq) -- 



2\l- Re / 

b 



\J (1 — v 2 )up\/l — sin 2 



2 ex P ( wr " \ ) } / 27TP 2 _ 4 sin(wr) 



2i/ - Re^= *• v A/ ' V Q 9 = i (18) 

V -~ %r 3 sin 2 0' / ~ 7* V ' 



7r / ~ / ; — - V ur 3 sin 2 0' / 
\ (1 — Vq)upJ1 — sin ury i 



Substituting this in Eq. Ijl5(l . introducing a new spherical system of coordinates, with the polar axis in z direction, 
and noting that 9' is the polar angle of the vector q gives 



7r— arcsin C/r 

Q> ™ a' 



, d6' sin f , „,,,n f , n sin9sin(ur) 

I(p,C,lj) = — -j- / = / rfgg 2 |/?M')l / dO- 



47r 4 r 7 / ~ 72 J ' J u J q 2 cos 2 Q + 2mU j 

arcsin C/r \/Slll 9' -0 



u = -q cos 6* + vV cos 2 + 2mw , 0(q,8') = b(qsin9' ,qcos9') , 
or, taking u as the integration variable, 

7T- arcsin C/r oo V q 2 +2?nu>+\q\ 

dO' sin 0' f, ,., /" , sin(ur) 



I ^M = -^r J 1 = = ]dqq\P{q,6Y J du^^ . (19) 

arcsinC/r J Sm 2 9' - — ^/q2+2mw-|<2l 

To further transform this integral, it is convenient to define a function T(q,9,uj) according to 



T(q,e) = ^2 / ,lw,r\ iKc.O'f. (20) 



Then integrating by parts, and taking into account that T(q,6,uj) — ► for q — * oo, brings Eq. (|19fl to the form 

7T — arcsin £/r +00 

me iut » f dO sin 6 f ^ r(g,0) 

arc: 

sin(w) 



I{P,C,u) = - 2 „ / ; / dq 



27r 2 r 7 /~ £2 J ^J q 2 + 2mw 

arcsin C/r y Sm " 2 ° 



sm(ur) 



u—\J q 2 -\-2muj — q 



(21) 



Finally, integrating / by parts once more, we find 



7T — arcsin £/r +oc 

~, , . e iut0 f d9sin9 f , 

I{PM = ^ I I e J * 

arcsin C/r W Sm ^ 2 ° 



x sin ^r-\/g 2 + 2r7iw^ |2rcos(gr) + -^-sin(gr)| . (22) 

This expression considerably simplifies in the practically important case of low ui and large r. Namely, if uj is such 
that 

D , , 

w < — = w , (23) 
mr 
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and also 

rD > 1 (24) 

(and therefore, oj <C D 2 /m), then the first term in / turns out to be exponentially small (~ e~ rD ) because of the 
oscillating product of trigonometric functions. Replacing sin 2 (gr) by its average value (1/2) in the rest of / gives 

7T — arcsin£/r +oc 7T— arcsin£/7" 

-. ^ . e luJt ° f d9 sinO f , dT e luJt ° f d6sm9T(0,( 



47r 2 r 2 tj J I ' £2~ J dq A-K 2 r 2 uo J I ' £2 



arcsinC/r \ Sin 9 5- arcsinC/r \/Sm 



. 7T-X +OO 

/(p, C ,a;) = --^- / - dds[n() { d qq \P(q,0)\ 2 , (25) 



x 



sm 2 — sin 2 \ q 



where x is the angle between the vector r and x, y plane, sin^ = C/ r - 

In the case of a spherically symmetric wave function, (3(q,6) — /3(g), integration over 9 in Eq. (|25() yields the 
expression derived in Ref. |lf)l | 



+00 

2 



I{pX^) = -^- J - J dqq\P{ q )\- (26) 


It has been assumed in the course of derivation of Eq. (|25|) that ui > 0, £ > 0. It is not difficult to verify that in the 
general case, I(p,(,u) = e^ u -\ u ^ to I(p, \(\, \cj\). 

Substituting the obtained expression into the defining equations Q14[l. ©, @, we thus obtain the following expres- 
sion for the low-frequency asymptotic of the correlation function 

T— IXI +OO 

C 00 (M, r » = e**-'>^^ J / **X,,»)P, (27) 

1 1 ^ y'sm 2 6'-sm J x 

all other components of the correlation function being suppressed by the factor \q\/m. 

In applications to microelectronics, u> varies from 10~ 6 Hz to 10 6 Hz, the relevant distances r are usually 10~ 5 cm 
to 10 _2 cm, D ~ h/d, where d « 10 -8 cm is the lattice spacing, and m is the effective electron mass, m w 10 _27 g, 
hence, lu w 10 10 Hz, so the conditions r » d, ui uj are always well-satisfied. 

In connection with the application of the obtained results to solids, it should be stressed that they refer to long- 
living free-evolving electron states. At the same time, because of collisions of electrons with phonons, impurities, 
and with each other, their evolution in a crystal usually cannot be considered free. It is important, however, that in 
view of smallness of the electron mass in comparison with the atomic masses, the electron collisions with phonons 
and impurities may often be considered elastic. Such collisions do not change the electron energy, and therefore 
they do not influence time evolution of the electron wave function. As to the electron-electron collisions, they do 
change the energy of electrons. However, the electron component in solids is practically always degenerate, and hence, 
only electrons with energies near the Fermi surface are actually scattered. Therefore, the above results concerning 
dispersion of the electromagnetic field fluctuations remain essentially the same despite the electron collisions, provided 
that they are applied to electrons far from the Fermi surface, and expressed in terms of the electron density matrix, 
rather than the wave function. Denoting the diagonal elements of the momentum space density matrix of electrons 
by p(g, 9), Eq. flTFI thus takes the form 

The power spectrum of the noise produced by uncorrelated electrons in a sample can be found by integrating Eq. 
with respect to Xq over the sample volume. If the time instants to are distributed uniformly (which is natural to 
expect), then the value of the total noise spectrum function, Cqq*, remains at the level of the individual contribution 
independently of the number of electrons in the sample, because of cancellation of the alternating phase factors 
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IV. THE NOISE INDUCED BY EXTERNAL ELECTRIC FIELD 



According to the property 3) of observed flicker noise, mentioned in the introduction, power spectrum of the noise 
induced by external electric field is proportional to E 2 , at least for sufficiently small values of the field strength E. If 
the density matrix g were analytic with respect to E, then the scalar function C}${E) — 0^(0) = AC^iE) would 
expand in even powers of E, and therefore, the leading term would be quadratic in the field strength, as required. 
However, g is not generally analytic in E, and therefore, the expansion of ACqq*(-E) might contain, e.g., a term 
proportional to \E\, in contradiction with the experiment. Thus, our primary concern below will be the weak field 
asymptotic of the power spectrum given by Eq. (|2"S)l . There is actually a simple and general reason why the right hand 
side of this equation should be quadratic in E despite possible non-analyticity of the density matrix g. This matrix 
is a functional of the equilibrium density matrix, go, and of the field strength. As was mentioned in the preceding 
section, the electron states responsible for the 1 jio behavior of the correlation function are those with energies far 
from the Fermi surface. For such states, the function g is a constant inversely proportional to the electron density, 
in particular, it is independent of parameters characterizing the Fermi surface, such as Fermi energy or momentum. 
Therefore, on dimensional grounds, the integral on the right of Eq. (|28l) should be proportional to the square of a 
characteristic momentum built from the field strength. The only such momentum is e|-E|r, where r is some time 
parameter characterizing electron kinetics, and hence, ACqq* ~ E 2 . This reasoning will be illustrated below using 
the simplest kinetic model with the relaxation-type collision term, which admits full theoretic investigation. 

In the presence of a constant homogeneous electric field, the model kinetic equation reads 

eEp- = -i^B. , (29) 
dq z t 

where z axis is chosen in the direction of E (E z = \E\ = E), r is the relaxation time, and e the electron energy. 
Since Eq. l|28|l was obtained for a free electron, we assume that the band structure in the given solid is parabolic, 
e = q 2 /2m*, where m* is the effective electron mass. According to Eq. <|II[) . the function g is normalized by 

(2^ (£) = 1 - 

If the function go satisfies this condition, then the normalized solution of the kinetic equation 129JI has the form 

oo 

' dte-V*>Qo(e'), e' = —( q 2 ± +e), <7o = \e\Er . (30) 



g = , „ 

qo J 2m 



<h 



Thus, in order to find the quantity Coo, we have to calculate the following integral 

7T — % +00 OO 

K= j ; dqq & - / d^<> g (e>) , S> = -^(q 2 sin 2 9 + £ 2 ) . 



X V Sm V~ Sm X ? cos0 

Changing the integration variables q, 9 — > £, r\ according to 



cos0=^===, q =C +V 



brings this integral to the form 



+oo , oo 

r r pCcosx/qo r , i 

K = JJdtdr ] 2<?o J dte-^g (e"), e" = —(rj 2 + C 2 sin 2 x + . 

— oo Ccosx 

Integrating by parts with respect to £ yields 

+ oo +oc oc 

— oo — oo ^ cos x 
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The first term in this expression represents the value of the quantity K in the absence of the external field. Therefore, 
performing a shift £ — > £ + £ cos x, and then £—>£ — £ cos \ in the remaining integral, we find the part induced by 
the electric field 



2 + °° 



— oo 

Going over to the polar coordinates in the (£, rj) plane gives finally 

oo +oo oo 

AK = 7T sin 2 x J rf£ £e~ 5/ <?° / de^ (e) = -tt sin 2 x J & 



£ 2 sin x 



£ 2 sin 2 x/2m* 



In the case of a degenerate electron system, the function Qo(s) is nearly constant up to Fermi energy where it falls 
off to zero. For such a function, AK ~ g 2 ,, up to exponentially small terms, provided that go is sufficiently small. 
Indeed, assuming go <C <?f, where gj? is the electron momentum at the Fermi surface, and neglecting terms of the 
order e~ qF l q ° , one can substitute go ( ^ x j by go(0) in the above integral to obtain 



oo 

AK = -7rsin 2 xeo(0) J d£ £ e -« / " = -gggo(0)7r sin 2 X ■ 





Inserting this into Eq. (|28() , we arrive at the following expression for the power spectrum of the induced noise 

327r 4 r 4 r'|o;| 

Thus, ACoo ~ E 2 , as was to be shown. It is important that the correction terms are of the order exp _<?F / 9 °, and 
hence, from the practical point of view, the condition go ^ If amounts to go < If- This implies that the upper limit 
on the field strengths for which Eq. i|31|) is valid is well above all experimentally relevant values. 

Let us compare this result with the power spectrum of quantum fluctuations in the absence of external field, 
obtained in Ref. |10|. which in the present notation has the form 



+oo 

327r 4 r 2 rVI 



C 00 (0) = e M*o-0 i dqqgo{£) _ 



o 

By the order of magnitude, 

Ag oo gogo(PJ „ n 2 n , m ll 

J dqqg 
o 

since Qo(0) ~ g^ 3 . Thus, in the model considered, the induced noise represents a relatively small correction. 

In practice, one is interested in fluctuations of the voltage, U, between two leads attached to a sample. Using the 
above results, it is not difficult to write down an expression for the voltage correlation function, C\j(x,x' ,u>). We 
have 

(U(x, x', t)) (U(x, x', t')) = (i (x, t) - A (x', t)) (i (x, t') - A (x', t')) 

= (A (x,t))(A (x,t')) + (A (x',t))(A (x',t')) 

- (A (x',t))(A Q ( Xl t')) - (A (x,t))(A (x',t')) , (32) 

where (• • • ) denotes averaging over the given in state. Fourier transforming and applying Eq. H31JI to each of the four 
terms in this expression, we obtain the power spectrum of the induced voltage fluctuation across the sample 

aC „ = -e^M(5! + £!_^_CL). (33) 
32-7T 4 \u>\ V T° r'° r'^r r^r I 
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This equation represents an individual contribution of an electron to the electric potential fluctuation. As was 
mentioned at the end of Sec. 1111 lH because of the oscillating exponent e luJt °, the magnitude of the total noise remains 
at the level of the individual contribution independently of the number of electrons. Therefore, summing up all 
contributions amounts simply to averaging over Xq : 

I a svtot I _ e 2 q%Q (0)G 

|A6(,| "l24#T' (34) 

where 

= /^o(5 + 5"S~^) (35) 
n 

is a dimensionless geometrical factor, and f2 is the sample volume. 

It is convenient to express the right hand side of Eq. I|34l) through the number of electrons in the sample, N, rather 
than the sample volume. Taking into account that for degenerate electrons, Qo(0) — 2/n, where n is the electron 
density (the factor 2 accounts for two spin states), substituting n = N/Q, qo = \e\Er, and restoring the ordinary units 
brings Eq. I|34|) to the form 

where / = |w|/27r. Finally, if the voltage leads are aligned in z direction, as is usually the case, Eq. can be 
rewritten as 

lA^toti uU 2 n e 4 r 2 G 

\ ACu I = W a = 52^W (37) 

where U is the average voltage across the sample, and I is the sample length in the direction of the electric field. 
In this form, it is similar to the well-known empirical Hooge law, except that N in Eq. Q37f) is the total number of 
electrons in the sample, rather than the number of charge carriers. We see that the analog of the Hooge constant, a, 
depends on physical properties of the sample material as well as on the sample geometry. 

Let us discuss the role of the sample geometry in somewhat more detail. Consider two geometrically similar samples, 
and let s be the ratio of their linear dimensions (see Fig.|3|). It turns out that such samples are characterized by the 
same value of the G-factor. Indeed, since the voltage across the sample is usually measured via two leads attached 
to its surface, the radius-vectors of the leads drawn from the center of similitude scale by the same factor s, and 
therefore, 

C 2 f ,3 (SZ - Z a ) 2 f q q (SZ - SZ ) 2 

r 4 r J \sx — Xo\ ^\sx ' — xq\ J \sx — sxo\\sx' — sxq\ 

o) 2 = [ d s x ? 

\x - x \*\x' - xq\ J °rV 

Oi Oi 

and likewise for the other terms in Eq. I|35|) . Thus, it follows from Eq. (|34|l that for geometrically similar samples, the 
noise level is inversely proportional to the sample volume. This is in agreement with the property 1) of flicker noise, 
mentioned in Sec. [I] It is also clear that the distribution of the noise magnitude around the value Cqq* is Gaussian, 
by virtue of the central limiting theorem. Thus, the quantum field fluctuations in a sample possess the property 2) 
as well. 

It is worth also to make the following comment concerning expression (l.'j.lti . The integrand in this formula involves 
the terms ( 2 /r 5 and £' 2 /r' 5 which give rise formally to a logarithmic divergence when the observation points approach 
the sample. In this connection, it should be recalled that the above calculations have been carried out under the 
condition rD 3> 1, hence, r, r' cannot be taken too small. Furthermore, one should remember that in any field 
measurement in a given point, one deals actually with the field averaged over a small but finite domain surrounding 
this point, i.e., the voltage lead in our case. Thus, for instance, the quantity Aq(x, t) appearing in the expression 1321) 
is to be substituted by 



^ = ?T / ^^o^'*)' 
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where J7l is the voltage lead volume. As a result of this substitution, the term ( /r , for instance, takes the form 



flf J J lai — £Eo| 4 |^ — a^ol 



S1 L 

Upon substituting into Eq. (|35[) . this term gives rise to a finite contribution even for intersecting f2, f^. 

V. CONCLUSIONS 

The main result of the present work is the general formula i|27|l describing the power spectrum of quantum elec- 
tromagnetic fluctuations produced by elementary particles in the presence of an external field. This formula shows 
that in the low- frequency limit, the power spectrum exhibits an inverse frequency dependence. Although the range of 
applicability of Eq. (|27|l . given by the conditions l|23|) . I|24|) . depends on the problem under consideration, it embraces 
virtually all experimentally relevant frequencies and distances. In particular, in application to microelectronics, the 
term "low-frequency limit" means that / <C 10 10 Hz, which covers well the whole measured band. To the best of the 
author knowledge, none of the other physical mechanisms of flicker noise, suggested so far, has been able to explain 
the observed plenum of the 1// law. 

We have applied the general formula to the calculation of the power spectrum of fluctuations produced by electrons 
in a sample in external electric field, assuming that the electron kinetics is described by a model equation with the 
relaxation-type collision term. This calculation shows that the power spectrum of induced fluctuations is proportional 
to the field strength squared for all practically relevant values of the electric field. We have also argued that this 
conclusion in fact holds true in the general case. Finally, we have established the exact dependence of the power 
spectrum on the sample geometry. We have shown that for geometrically similar samples the noise power spectrum 
is inversely proportional to the sample volume, while for samples of the same volume dependence on the sample 
geometry is described by the dimensionless G-factor given by Eq. i|35[l . 

Qualitatively, the established properties of quantum electromagnetic fluctuations match perfectly with the ex- 
perimentally observed properties of l/J-noise, and suggest that these fluctuations can be considered as one of the 
underlying mechanisms of flicker noise. As to the quantitative side, the estimates of Ref. 01 show that the noise 
level predicted by Eq. I|27|l in the absence of external field is in a reasonable agreement with experimental data. At 
the same time, according to the results of Sec. IIVI the change of the noise level in external electric field is relatively 
small, which disagrees with observations. Of course, it is difficult to expect that predictions of the simple model 
employed in Sec. IIVI will be quantitatively correct. Although the found disagreement may be the result of the model 
oversimplification, it raises the question of possible alternative mechanisms of the noise amplification by external field. 
Recall that the expression l|34(l for the power spectrum of fluctuations produced by a sample was derived assuming 
uniform distribution of the time instants to- Therefore, correlation between irj's is a possible source of the noise am- 
plification. Since the number of electrons in the sample is large, already a relatively small correlation in the values 
of to would result in a noticeable increase of the noise level. Thus, the question concerning the relative role of the 
quantum electromagnetic fluctuations in explaining the observed 1// noise requires further investigation. 
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FIG. 1: Feynman diagrams representing the disconnected part of correlation function [the second term in Eq. @]. Wavy 
lines denote the photon propagators, solid lines the massive particle, q and p are the particle 4-momentum and 4-momentum 
transfer, respectively. 




(a) (b) (c) 

FIG. 2: Feynman diagrams representing the connected part of correlation function [the first term in Eq. P)]. 




FIG. 3: Voltage measurement in geometrically similar samples (s = 2). 



